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Abstract. We compute the charm drag and diffusion coefficients in a hot pion gas, such as is formed 
in a Heavy Ion Collision after the system cools sufficiently to transit into the hadron phase. We fully 
exploit Heavy Quark Effective Theory (with both D and D* mesons as elementary degrees of freedom 
during the collision) and Chiral Perturbation Theory, and employ standard unitarization to reach higher 
temperatures. We find that a certain friction and shear diffusion coefficients are almost p 2 -independent at 
fixed temperature which simplifies phenomenological analysis. 

At the higher end of reliability of our calculation, T ~ 150 MeV, we report a charm relaxation length 
A c ~ 40 fm, in agreement with the model estimate of He, Fries and Rapp. 

The momentum of a 1-GeV charm quark decreases about 50 MeV per Fermi when crossing the hadron 
phase. 

PACS. 14.65.Dw Charmed quarks - 25.75. Ag Global features in relativistic heavy ion collisions -51.20.+d 
Viscosity, diffusion, and thermal conductivity - 12.39.Fe Chiral Lagrangians - 12.39.Hg Heavy quark 
effective theory 



1 Introduction 

Heavy Ion Collisions provide a thriving branch of nuclear 
and particle physics. Thanks to technological advances in 
the last three decades, measurements that once looked too 
challenging can now be performed. One of these is the 
reconstruction of charmed and bottomed mesons flowing 
out of the nuclear debris, that modern vertex detectors, 
together with good particle reconstruction and the abil- 
ity to automatically treat very large data samples have 
brought to the realm of measurability. 

Heavy-flavored hadrons are interesting because the 
hadron medium is not hot enough to excite charm pairs. 
They are produced by hard gluons in the initial stages 
of the collision and their spectra will carry a memory of 
it, unlike pions and kaons that can be produced in the 
thermal medium at later stages, and thus show a spectrum 
close to black-body without much information from the 
initial configuration of fields. 

However, charmed and bottomed mesons do interact 
with the hadron gas after the crossover from the high- 
energy phase (that, although now known to be strongly 
coupled, we will continue naming "quark-gluon plasma" 
as is customary) . The corrections to their properties due to 
this cooler medium requires their scattering cross-section 
with the medium pions and other particles. Given the scat- 
tering amplitudes one can proceed to kinetic simulations 



following individual particles, or employ kinetic theory to 
compute transport coefficients that can be input to bulk 
hydrodynamic simulations. 

In this article we will be concerned with charmed mesons, 
of more immediate interest, although the theory devel- 
oped can immediately be applied to bottomed mesons too, 
which we will leave for a future application. 

The scattering amplitudes or cross sections for heavy 
mesons cannot be directly accessed by experiment (since 
the short life of these mesons makes impossible to fo- 
cuse beams of them on a target) so their knowledge re- 
quires theory constraints. In the past [1] cross-sections 
were only guessed on the basis of constituent quark count- 
ing. Since a pp is about 40 mbarn and ct^n about 2 mbarn, 
this counting leads to a charm quark-light quark scatter- 
ing cross-section which is a cq ~ 0.3 mbarn, much smaller 
than a qq ~ 4 mbarn, leading to ao n ~ 9 mbarn. 

As we will see, this old reasoning is not too much off 
the mark, but the cross-sections can now be accessed with 
more reliable theoretical methods 0EJE1EJIMI1IE] j combin- 
ing Chiral Perturbation Theory, Heavy Quark Effective 
Theory, and Unitarity. 

Given the renewed experimental interest, it appears 
that several theoretical groups have simultaneously been 
attempting to extract the transport coefficients from the 
increased understanding of hadron-hadron interactions. 
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1.1 Current theoretical understanding and setup 

The work of Laine [9] employs canonical perturbation the- 
ory in HQET and ChPT and thus focuses on the lowest 
possible temperatures. Two simultaneous papers of He, 
Fries and Rapp [TO] and of Ghosh et al. [H] have at- 
tempted to reach higher temperatures, close to the cross- 
over to the quark and gluon plasma, by including further 
species of particles (K and r\ mesons or nucleons). While 
the second combines the perturbative approach of Laine 
with Born exchange terms, the first relies back on phe- 
nomenological estimates of the cross-sections. 

We feel that there is still room for our contribution. 
There are serious disagreements among the three works 
cited. The extension of Ghosh et al. to higher energies does 
not include unitarity as a guiding principle, thus likely 
overestimating the cross-section since the polynomial per- 
turbative expansion grows very fast with s. 

By performing a state of the art computation of the 
pion-charmed hadron interactions, extending the work of 
Laine and Ghosh et al. by providing both the canonical 
HQET+ChPT perturbation theory and unitarity, and ty- 
ing the unknown parameters to experimental Dq and D\ 
resonances, we believe we have an interaction that is both 
solidly grounded in theory, and phcnomenologically ac- 
ceptable, drawing from the best features of the extant 
works. 

As the charm transport coefficients are concerned, we 
will consider the drag or friction force F (variously de- 
noted 7, r\ or A in the literature), and the two -To and F\ 
momentum-space diffusion coefficients. Other works have 
considered only isotropic drag and diffusion, in which case 
there is only one diffusion coefficient also denoted as K or 
B . We do not make this hypothesis and provide both co- 
efficients corresponding to parallel and shear momentum 
transfers. Finally, in the p — > limit, we make contact 
with the traditional kinetic theory and compute the space 
diffusion coefficient D x (again, sometimes denoted D s in 
the literature, but we avoid this notation to prevent con- 
fusion with the meson of equal name). We find important 
to lift the hypothesis of isotropy because of the interesting 
elliptic flow observable. 

Finally we make an additional contribution, in the phi- 
losophy of fully exploiting Heavy Quark Effective Theory 
(in addition to ChPT) as a starting point. In the Heavy 
Quark limit, the D and D* mesons are degenerate and 
there are four (one spin-zero and three spin-one) prop- 
agating modes for the charm quark in the pion medium. 
This has been missed by all existing approaches, that only 
attend to D-propagation in-medium (Laine however con- 
siders both B and B* in the bottom sector, where the B* 
meson is stable under strong decays). 

In the physical world the D* meson is unstable and de- 
cays to Dir. However it does so with a small width (given 
its closeness to threshold) and thus, for the duration of the 
hadron gas, expected to be of the order of 5 — 10 fm/c at 
most, it propagates as a stable mode. To settle this point 
let us realize that, for a particle to decay within 5 fm of 
its production point, its width has to be of order 40 MeV, 
which, even after accounting for in-medium modifications 



(see Fig. 1 of [TO]), is only reached for D* mesons at tem- 
peratures of order the phase transition T ~ 180 MeV, so 
that for the entire life of the hadron gas both D and D* 
mesons need to be taken as elementary degrees of freedom. 

We of course include the DD*ir interaction vertex in 
the effective Lagrangian. However we will present compu- 
tations in which the D* is thus included as an elementary 
particle (but also others without it for ease of compar- 
ison with the recent computations). We will gcnerically 
speak of the passage of the charm quark through the pion 
medium, whether hadronized in a D or a D* meson. In the 
heavy quark limit, the heavy quark is little affected by the 
specific nature of the light degrees of freedom hadronizing 
around it. 

Given these theory improvements, and the fact that 
the other groups have not found very large effects from 
including strangeness or nucleons as explicit degrees of 
freedom in the hadron gas, we will comptent ourselves 
with examining the contribution of pions. A priori one can 
expect pions to provide the bulk of the charm-medium in- 
teraction, by their large multiplicity (typically one particle 
of any other species for every ten pions). 

We employ the Fokker-Planck formalism for a heavy 
Brownian particle subject to the bombardment of the light 
pions in the medium. Our approximations will be sensible 
as long as the momentum of the heavy particle remains 
smaller than its mass in natural units, so that p > 2 GeV 
is not accessible by our computation (although we show 
plots at higher momentum for ease of comparison with fu- 
ture investigations addressing hard heavy flavors). Pairs 
of heavy quarks rapidly drift apart and, by the time of 
the transition to the hadron phase, they are at least three 
Fermi away from each other and never rescatter in it (un- 
less initially in a bound charmonium state). Since they 
are very scarce, we neglect the interactions between charm 
pairs formed in different points of the collision. 



1.2 Experimental motivation 

One could conceive hydrodynamic calculations of the quark- 
gluon plasma that would result in fits of F, Fq and A to 
experimental data under certain assumptions on the ini- 
tial distribution of heavy quarks. The information gained 
would be very valuable to understand how strongly that 
plasma is coupled, and perhaps restrict the possible initial 
state configurations. 

However, the extraction of the coefficients is blurred 
by the hadron phase in the final state, as the system must 
cool before total freeze out, and charm quarks propagat- 
ing through the resulting hadron medium will also suffer 
drag and diffusion. It is like trying to deduce the dispersive 
properties of a glass with a beam of light going through 
an additional lens: both have to be simultaneously under- 
stood. 

Existing data on nuclear suppression factors and ellip- 
tic flow (see Sec. [6] below) have already been compared 
with standing calculations within the asymptotic quark- 
gluon plasma phase [12) and, perhaps more successfully, 
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with a mixed approach that includes resonances surviving 
into the plasma phase |13j . 

Another observable that is being addressed in the lit- 
erature is the transverse momentum spectrum of the D 
mesons, that should be a rough thermometer of the phase 
transition fj], provided that the effect of the final stage 
hadron phase does not blur all information out (it doesn't, 
as we will show in this article). 

1.3 The £>-meson spectrum 

A charm quark propagating in the low-temperature medium 
below the deconfinement phase transition must do so con- 
fined in a hadron. In central heavy-ion collisions the baryon 
number is very small and can be neglected. Therefore one 
expects the charm quark to form a Z)-meson or an ex- 
citation thereof. Let us briefly recall what experimental 
knowledge there is about the D-spectrum. 

The ground state D-mcson is as usual in meson spec- 
troscopy a pseudoscalar J p = 0~ with four charge states 
+, — , 0, (identified in the quark model as cd, dc, cu and 
uc respectively in a relative s-wave with spins antiparal- 
lel). Since we neglect isospin-breaking terms, we can av- 
erage the masses over this quartet to obtain Mo ~ 1867 
MeV. 

This meson cannot decay by any strong process and we 
will take it to be absolutely stable. 

The first excitation is the vector 1 _ D* meson whose 
mass average is Md* = 2008.5 MeV. In the Heavy Quark 
Limit this meson should degenerate with the D, (and in 
fact this is seen by glancing higher to the i?-meson whose 
splitting to the B* is much smaller). This mass is barely 
above Dtt threshold, so there is only this one strong decay 
channel, and it is very suppressed. 

The width of the charged D* is estimated at 1 MeV, 
and that for the neutral partners has not been measured 
but is consistent with r < 2 MeV. This means that a 
D* has a mean lifetime in vacuum of order 100-200 fm. 
Since the typical freeze-out time of a heavy ion collision 
is about 20 fm it is not a bad first approximation to take 
the D* meson as also stable during the fireball's lifetime: 
there is ample room even in medium since the decay time 
is an order of magnitude larger than the freeze-out time. 
As stated above, in-medium corrections do not alter the 
picture. Thus, most D* mesons decay after collisions have 
ceased. This approximation can be corrected if wished by 
taking into account the in-medium inelastic process D* — > 
Dn with Bose enhancement for the final-state pion. 

In agreement with quark model expectations, the next- 
higher excitations of the D system seem to be a triplet 
and a singlet of positive parity, with spins + , 1 + , 2 + and 
1 + respectively, corresponding to 2S+1 Lj = 3 Pj and 1 P\. 
The two mesons with spin 1 and positive parity must mix, 
and they do so in an interesting manner: the one with low- 
est mass, £>i(2420) becomes narrow and hence decoupled 
from the natural s-wave decay channel D*ir, whereas the 
higher member £>i(2430) is very broad and seen in that 
configuration. The situation can be seen in Fig. [1] and in 
Table Q] 




D 2 (2460) 



71 

/ 

D*(2008) 



D (1867) 



Fig. 1. The currently known low- lying D- meson system. The 
negative parity state D and D* are represented as the blue 
lines. The four positive parity states have the mass measure- 
ment spread throughout the red boxes, while the hollow black 
boxes represent current estimates of their width, s-wave pion 
decays are depicted. 



Table 1. Charged-average masses and experimental esti- 
mates Q3] for the strong widths of the D-meson resonances. 
Units are MeV. Errors not quoted are about 1 MeV or less. 



Meson 


J p 


M (MeV) 


r (MeV) 


D 


0" 


1867 




D* 


1" 


2008 


1 


D 


0+ 


2360(40) 


270(50) 


D 1 


1+ 


2422 


22(5) 


D 1 


1+ 


2427(40) 


380(150) 


D 2 


2+ 


2460 


30 



The remaining low- lying resonance, the Z?2, is again 
narrow. Since its mass at 2460 MeV is 600 MeV above the 
ground-state D meson, and it is quite decoupled due to 
its moderate width of about 40 MeV, we do not expect 
this (nor the Z?i(2420) to play an important role at small 
temperatures. 

Thus a sensible approach to charm propagation in a 
heavy-ion collision after the phase transition to a hadron 
gas has occurred, is to take the D and D* mesons as ab- 
solutely stable degrees of freedom for the c-quark, that in 
collision with the in-medium pions they rescattcr into the 
resonances D and Di(2430). 

The experimental knowledge of the resonances Dq and 
Di sufficiently constrains the low-energy effective Lagrangian 
density for cir scattering so that we are in possession of a 
good approximation to the cross section. 
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2 The Fokker-Planck equation 

2.1 Derivation 

The momentum-space distribution of charm quarks with 
momentum p, / c (p), is not in equilibrium when the hadron 
phase of a heavy-ion collision forms, and must relax via a 
Boltzmann equation. 



dt 



C[/ C (P)] 



(1) 



The right hand side is called the collision operator, be- 
cause it describes kinetic collisions of the charmed parti- 
cles. The left hand side, in the absence of external forces, 
is the advective derivative 



dfdp) 
dt 



v • V x / C (p) 



dfc(p) 



dt 



(2) 



coll 



The density of D and D* mesons being very small, 
we can neglect collisions between D mesons themselves 
and concentrate only on the interaction of these charmed 
mesons with the pion bath, assumed in thermal equilib- 
rium. 

The bath's distribution function / ff (q) is hence the 
Bose-Einstein function. Moreover, the gas is assumed ho- 
mogeneous and the distribution does not depend on x. For 
this reason one can average the Boltzmann equation over 
the collision volume and understand the 1-particle distri- 
bution function for the charmed mesons as the average 



/c(*,p) = y J dx / c (t,x,p) 



The averaged Boltzmann equation becomes then 



df c (t,p) 
dt 



df c (t,p) 
dt 



(3) 



(4) 



coll 



Charmed mesons may enter and exit the momentum 
element dp around p by collisions with the pion bath, so 



the collision term has two parts associated with gains and 
losses. 

Gains in the momentum distribution around p are 
proportional to the probability density around (p + k) 
times the probability of transferring momentum k from 
the charmed meson to the bath. It is therefore convenient 
to define a collision rate w(p, k) for a charmed meson with 
initial and final momenta p, p — k. 

Conversely, losses are proportional to the distribution 
function around p times the probability of transferring 
momentum k to the pion bath. 

In principle, the Boltzmann equation should be treated 
as a quantum Boltzmann- Uchling-Uhlcnbeck equation tak- 
ing into account Bose enhancement effect in the final state, 
with factors (1 + f c ) that encode the increased probability 
of a charmed meson scattering into an already occupied 
state, 



df c (t,p) 
dt 



(5) 



J dk{/ c (t, p + k)w(p + k, k) [1 + f c (t, P)] 



-/c(t,pMp,k)[l + / c (f > p-k)]} . 

However, as the number of c-quarks is very small, we 
can approximate 1 + f c (t, p) ~ 1 inside the collision oper- 
ator in practice. This approximation however is probably 
not valid for the pion distribution function and we keep the 
(l+f„) factor in Eq. ([7]) below. As the charmed mesons are 
concerned, a classical Boltzmann equation should however 
be very accurate, 



df c {t,p) 
dt 



= Jdk [/ c (t,p + k)ffl(p + k,k)-/ c (t 1 pHp,k)] 

(6) 

In turn the collision rate can be spelled out in terms of the 
Lorentz invariant charm quark-pion scattering amplitude, 



w(P,k) = fjrr 



dq 



IM) [l + / ff (q + k)] 



1 1 



2E- 2E° 2E- +k 2E. 



p—k 



-(2n) 4 S(E c p 



El 



E 



p—k 



E q+k) 



\M vc (a,t, X )\ 2 
(7) 



(<7tt = 3 is the pion isospin degeneracy, and x denotes 
the possible spin degrees of freedom, active if the c quark 
finds itself inside a, D* meson). The scattering amplitude 
A4 is normalized according to standard covariant con- 
vention [H]. Note that Eq. (7) of PU] differs by the 
Bosc-enhancement factor (1 + f T ) for the pion exiting 
the collision. We believe that in the temperature range 
of ~ T ~ 150 McV that we (and those authors) treat, 
this enhancement should not be neglected. 

The Boltzmann equation in this case reduces to a much 
simpler Fokker-Planck equation because the mass of the 



D and D* mesons carrying the c-quark is much greater 
than the mass of the pions and the temperature of the 
heat bath. Then, the scale of momentum for which there 
is a significant change of f c {p) with the momentum of 
the D meson |p| is greater than the typical transfcred 
momentum |k|, that is of the order of T: 

|p| fc > |k| ~ T ~ 150 McV . (8) 

Because of this separation of scales, it is natural to ex- 
pand the collision rate inside the collision operator respect 



Abreu, Cabrera, Llanes-Estrada and Torres-Rincon: Charm diffusion in a pion gas.. 



5 



to its first argument p + k. 



wf = w(p + k,k) / c (i,p + k) = (9) 
d 2 



(Wf) • 



with i,j = 1,2,3. The collision integral reads, with this 
substitution, 



\df c (t, P y 


= f dk 


dt 


coll J 



_d_ 

dpi 



if- 



2 ktkj dpidpj 



This suggests defining two auxiliary functions, 
r ij(p) = \ J dk w (P> k ) k ikj , 



(10) 

(11) 
(12) 



whose classical interpretation is that of a drag force act- 
ing on the charmed particle, and the autocorrelation of a 
random, Brownian force, as will be shown below in Ap- 
pendix [Bj 

Eq. (O reduces to the Fokker-Planck equation 



a/ c (*,p) 
at 



d_ 

dpi 



Fi(p)/c(* ) p) + ^[itf( P )/ c (t,p)] 



(13) 

where meanwhile we can see that Fi behaves as a friction 
term representing the average momentum change of the D 
meson and Zy acts as a diffusion coefficient in momentum 
space, as it forces a broadening of the average momentum 
distribution of the D meson. This interpretation also falls- 
off from the one-dimensional solution that we leave for 
Appendix [X] 

We will not find necesary to solve the three-dimensional 
Fokker-Planck equation for f c (p) in full, but only to cal- 
culate the coefficients Fi and F^ that already encode the 
physics of charm drag and diffusion. 



2.2 Fi and F^ Coefficients 

In the ideal case where the pion gas is homogeneous and 
isotropic, and because the coefficients Fi and only de- 
pend on p l , they can be expressed as a function of three 
scalar functions by means of 



Fi(p) = F(p 2 ) Pi , 

r ij { P ) = r ( P 2 )A ij + r 1 (p , 



P 2 



where 



PiPj 

r,2 



satisfies the handy identity AijA 



V 

ij = 2. 



(14) 



(15) 



We choose the momenta of the elastic collision between 
a charmed meson D or D* and a pion as 



D(p) + 7r(q) -> D(p - k) + 7r(q + k). 



(16) 



The three scalar coefficients in Eq. (TT41 are then simple 
integrals over the interaction rate 



F{p 2 



P 1 F, 



dk w(p,k) -Ar- 



(17) 



r Q ( P 2 ) = = I J dkw( P ,k) 

r 1 (p 2 ) = ^r i j = ydkw(p,k) 



k 2 

p2 



(hp^f 



where the dynamics is fed-in by the scattering matrix el- 
ements I-MttcI- The choice of kinematic integration vari- 
ables and the reduction of these integrals is detailed in 
Appendix [D] 

We also remind the reader in Appendix fB] how the 
interpretation of the- friction coefficient times the quark 
momentum F p is that of an energy loss per unit length 
upon propagation of the charm quark in the plasma, and 
how the loss of momentum per unit length is simply F E 
in terms of energy and momentum of the charmed particle. 

After we numerically control the cross-section and scat- 
tering amplitude M for the charm quark in the pion medium, 
we evaluate the three transport coefficients and give the 
results in subsection 15.21 below. We quote there two differ- 
ent approximations. One in which the D* is neglected as a 
propagating degree of freedom (akin to what can be found 
so far in the literature) , and one in which the c-quark can 
travel also as a D* meson (with slightly modified interac- 
tion and kinematics). 



3 Effective Lagrangian for D, D* and n with 
ChPT and HQET 

Now we construct the chiral Lagrangian density that de- 
scribes the interactions between the spin-0 and spin-1 D- 
mesons and pseudoscalar Goldstone bosons. The leading 
order (LO) chiral Lagrangian £W is given by [SHHfS], 



+m 2 n D*^D*} 



{D*^u n D ] - Du^D*}) 



2m u 



(18) 



where D = (D°,D+,D+) and D* = (D*°, D*+ , £>*+) M 
are the SU(3) anti-triplets of spin-zero and spin-one Ja- 
mesons with the chiral limit mass m^, respectively. We 
have also used the quantities 



where 



m ( m 

u = VU = cxp I - 



\V2F 



(19) 
(20) 
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is the unitary matrix incorporating the pscudoscalar Gold- 
stone bosons, 

(fr' + jm * + K+ \ 

F in Eq. ([20]) is the Goldstone boson decay constant in 
chiral limit. 

The NLO chiral Lagrangian L^ 1 reads 

C {2) = —hoDD^ < X +> +h 1 D X +D^ + h 2 DD^ < u^u p > 
+h 3 Du f "u fl D^ + h 4 \J ll D\J u D' [ < u^u v > 
+h 5 V p D{u^u"}V u D^ + h D^D*J < X +> 
-hD^x+Dp - h 2 D*»D*J < u v u v > 
-h 3 D*^u v u u D*J - hiV^D** V V B$ < u^u v > 
-hV p D* a {u^u v }V v D*J . (22) 

where 

X+ = ^X"^ + U X U , (23) 
with x = diag(mj, toJ, 2/71^ — m 2 ) being the mass ma- 
trix. The twelve parameters hi, hi(i = 0, 5) are the low- 
energy constants (LECs), to be determined. However, we 
can make use of some constraints to reduce the set of free 
LECs. First, it should be noticed that in the limit of large 
number of colors (N c ) of QCD [16], single- flavor trace in- 
teractions are dominant. So, we fix ho = h 2 = h± = ho = 
h 2 = hi = henceforth. Besides, by imposing the heavy- 
quark symmetry (as will become clear in subsection 14.11 
below), it follows that hi ~ hi. 



In the following, the lowest order of the perturbative 
expansion of the quantities i"^, and \+ m Eqs. (|18|) 
and (|22p is considered to construct the scattering matrix 
of the interactions between the charmed mesons and the 
pscudoscalar Goldstone bosons. 

4 Scattering matrix for the c quark in the 
pion gas 

From the Lagrangian in Eq. (|18[l we are able to obtain the 
scattering amplitudes V for D,D*<P — > D,D*<P processes. 
In Fig. [2] we show the tree-level diagrams constructed from 
the LO and NLO interactions. These include both contact 
interactions and Born exchanges. The different scattering 
channels are labeled as V a through V4, where the sub- 
scripts refer to the scattering channels as follows 

(a) : Dcf) -t Dcf), 

(b) : D*(j) -> D4> 

(c) : D(j) -> D*(t) 

(d) : D*<f> -> D*cj> . (24) 

Notice that, because of the scarcity of strange quarks 
in the heavy ion collison debris (kaon multiplicity is 10% 
of typical pion multiplicity) we are interested only in chan- 
nels involving the scattering between charmed mesons and 
pions with total strangeness equal to zero. So then we can 
simplify <f> it and write down the relevant amplitudes 
as 



Va = ~ + ~^T JLhl + -pT h 3(P2 - Pi) + ^T^5 [(Pi ■P2)(P3 - Pi) + (Pi 'P4)(P2 ' Pa)] 



'-prP 2 [CiD^ip! +p 2 ) + C 5 D^(p 2 - p 3 )]P4 



v b = 



ig z 



.2 



C 4 p£ (2pf + pf) P 4 P D»r( Pl + P2 ) + C 5 pt (pf -pi -pf) P2 P D v "ip 2 -ps)] e afj ^{ Pl ) 



Vr = 



m^F 2 



C 4P 2 (fx+P^+pl) P2 P D^( Pl +p 2 ) + C 5 p% U - 2pl) P4p D»?(p 2 -p 3 )] e a ^ v e*»(p 3 ) , 



V d 



~ |^( s + 2C ^r ?r hl + ^rMP2 -Pi) + ^T h 5 [(Pi -P2)(Pi -Pi) + (Pi -P4)(P2 -P3)]| e A1 (pi)e*(p 3 ) 

+ ^f- [CiDfa +P2)+C 5 £>(P2-P3)]P 2 %(P1K<(P3) 



Iff 



3m 2 D F 2 



C 6 p a 2 (2p£ + pf) pi (pi +pl+ pf ) D^( Pl + pa) + C lP a 2 (pf - 2pf ) p£ (p£ - p£ - p?) D^(p 2 - p 3 ) 



(25) 



where Cj (i = 0, 7) arc the coefficients of the scatter- 
ing amplitudes for D?t, D*tt channels with total isospin J, 
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Table 2. Coefficients of the scattering amplitudes for the 
Dn,D*n channels with total isospin / in Eq. (|25|) . 



Constants 


'=* 


1= * 
1 2 


Co 


-2 


1 


Ci 


-1 


-1 


c 2 


1 


1 


c 3 


1 


1 


c 4 


3 





c 5 


i 

3 


2 
3 


c 6 


3 





cv 


i 

3 


2 
3 



Fig. 2. Tree- level diagrams relevant to the scattering ampli- 
tudes in Eq. (|25[) . The solid, double and dashed lines represent 
the D-mesons, D*-mesons and Goldstone bosons, respectively. 

done in Tabled and D(p), D^ v {p) are the propagators of 
D and Z?*-mesons, respectively, 



D(p) 



— m 
—i 



P 



m 



p»p u 



D 



(26) 



As the two particles in all amplitudes are distinguishable, 
there is no t-channel type contribution (as e.g. in Compton 
scattering) with our relevant fields (open charm mesons 
and pions), and only s and it-channel interactions appear. 
Between a D and a 7r one could exchange additional, closed 
flavor resonances in the ^-channel, but a quick examination 
makes clear that these contributions are totally negligible. 
For example, /o exchange, while having strong coupling to 
two pions, has negligible coupling to two D mesons, so one 
of the vertices makes the amplitude very small. Similarly, 
J/tjj i-channcl exchange is suppresed because of the small 
two-pion coupling of the very narrow state (and similar 
for other, closed flavor resonances). It doesn't make sense 
to include these resonances while neglecting higher order 
chiral and heavy quark corrections to the Dtt Lagrangian 
with the basic fields. 

Finally e M (p) is the polarization vector of the vector 
D*-meson. If we were to write the polarization indices 
explicitly, e"(p) = e£(p), V b = V bx> V c = V cX , V d = W, 
while V a remains a scalar as no vector mesons appear. 

The amplitudes V4 and V c must be related by time 
reversal, since they encode D*tt — > Dtt and Dtt — > D*tt 
respectively. Indeed, if one exchanges p\ by ps and pi by 
Pi, and employs energy-momentum conservation pi +p2 = 
P3 + f>4) they map onto each other as V> — > V c , V c — > Vb- 



4.1 Heavy Quark Symmetry 

Manifesting the spin of the vector meson, but ignoring the 
isospin index, Heavy Quark Symmetry should manifest 
itself [9] by the transformation 



SD = -a ■ D* 
SB* = aD + axT>* . 



(27) 
(28) 



The Lagrangian density in Eq. (f2"2"|) has been constructed 
manifestly maintaining chiral symmetry, that is then bro- 
ken only carefully in perturbation theory upon expanding 
in fields and derivatives to construct the LO and NLO 
chiral amplitudes. 

However since the charmed quarks are heavy fields, 
one should recover the Heavy Quark Symmetry by taking 
mjj — > oo. 

Referring to our amplitude in Eq. (|23|) . both V> and 
V c are of order l/mo- To see it, one needs to write the 
denominator of the propagator as 



(p?- m \ 



((p° - m D ) x 2m D )~ 



and notice that in the numerator at most one of the mo- 
menta can take the value p° ~ mjj because of the anti- 
symmetric Levi-Civita tensor: the other three four- vectors 
have to take spatial values to avoid its vanishing. There- 
fore the term within brackets is of 0(1) in the 1/tod 
counting and the explicit factor of 1/mjj in front of the 
bracket suppresses the term. 

Thus, the spin-changing amplitudes D*tt — > Dtt and 
Dtt — > D*tt are of 0(1/tod) and vanish in the heavy quark 
limit as expected, a collision with a pion cannot change 
the heavy-quark spin that decouples. 

Turning to the elastic D*tt — > D*tt amplitude, we no- 
tice that the last bracket of Vd (carrying terms propor- 
tional to Cq and CV) is also NLO in the heavy quark 
counting. We see that because of the two Levi-Civita e 
tensors, only one of the pair of a, /3 indices and only one 
of the pair p, a can take the value 0. Because of the ex- 
plicit propagator of order I/tud, the bracket is at most of 
order m^, and the l/m 2 D factor in front of it suppresses 
it. 

The remaining part of the Vd amplitude is then equal 
to V a if we impose hi = hi as we have, yielding the ex- 
pression of Heavy Quark Symmetry 

V a (DTT -> Dtt) = -V d (D*n -> D*tt) , 

(the polarization of the D* will be handed shortly). That 
is, in the infinitely heavy quark limit, the charmed mesons 
propagate in four states (the D and the three polarizations 
of the vector D*) that do not mix with each other and have 
equal scattering rates with the pion gas. 

In this limit, the masses mo = m* D and also the dy- 
namical resonances accessible in the scattering have equal 
mass 7Ti£> = mjj 1 and width _Td = rjj 1 . 
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Further, the Born exchange piece (terms proportional 
to C4 and C5) in both Dir and D*ir elastic amplitudes is 
subleading in HQET. To demonstrate it, we expand the 
intermediate meson propagators 



iiV(pi+Pa) = (29) 



{Pi +P2) 2 - m\ 1 



(P?) 



2p\p\ 



„0 



1 



suppressed by m~j^ as claimed. 

The sum over polarizations also simplifies in the heavy 
quark limit. The vector-meson polarization basis then be- 
comes a conventional spacelike spin-1 basis tied to a fixed 
reference frame, in the Cartesian basis simply 



^(p, A)~e"(A) = (0,e A ) 



(30) 



satisfying the closure relation (/j, = and v = no more 
contributing) 



^e»(A)e-'"(A)=5 i 



A=l 



and an orthogonality relation 



i (X 1 )e»(X 3 ) = -6 XlX3 



(31) 



(32) 



with the minus sign from the spatial part of the metric. 
This sign cancels the explicit sign in front of the brace 
of the first line of Vd in equation (|2"5)) . Thus, the final 
amplitude for scattering off a heavy quark in the pion 
gas, to next to leading order in the chiral expansion and 
leading order in the heavy quark expansion, irrespective 
of whether the heavy quark is in a D or a D* meson, is 
given by 



2CW 2C 2 
- "y t p2 hi + j^-h 3 {p 2 ■ pa) + (33) 

2C 3 

-p2 h 5 [(Pi -P2)(P3 - Pi) + (Pi -Pi){P2 -Pz)] ■ 



The I/mo pieces included in our amplitude Eq. (|25[) 
are of course not the complete amplitude of NLO heavy- 
quark counting. It includes only those that are simulta- 
neously NLO in the chiral counting; we are dealing with 
a double series expansion of the total amplitude. Due to 
those corrections, and also for our allowing the physical 
D, D* masses to be slightly different, the properties of the 
Dq and D\ are not precisely the same. However they are 
close enough for most purposes. 



and, in our case, Born terms with an intermediate propa- 
gator due to the DD*ir coupling) and by their very nature 
are unable to describe excited elastic resonances (in our 
case, D and D±), 

The key to understanding this limitation is to note that, 
at fixed order, ChPT violates unitarity as momentum is 
increased. Therefore several strategies have been adopted 
to bypass the shortcoming, such as the N/D method, the 
Inverse Amplitude Method, or the K-matrix method. 

We pursue the simplest partial-wave unitarization by 
employing on-shell factorization |17j which is a nice fea- 
ture of polynomial expansions and leads to algebraic for- 
mulae for the unitarized partial wave amplitudes, capable 
of reproducing resonances. Our conventions for the expan- 
sion of the perturbative V a and unitarized T a amplitudes 
in terms of Legendrc polinomials arc 



VI = 



dx Pi(x)V a (s,x) 



dx Pi(x)T a (s,x) 



(34) 



(35) 



where x = cos 9 and Pq (x) = 1 and a is a channel index. 

We proceed by projecting the perturbative amplitude 
into the s-wave, that dominates at low energies because 
of the k 2l+1 suppression of higher waves, and is resonant 
at the Do (for Dir scattering) and D\ (for D*ir scatter- 
ing), thus dominating the entire amplitude at moderate 
heavy-quark velocities (at higher velocities, boosting to 
the moving center of mass frame kinematically induces 
higher waves) . Thus the perturbative amplitude is substi- 
tuted for 



dxV a {s,t(x),u(s,t(x))) P (x) . (36) 



The unitarized scalar amplitudes T a decouple in leading 
order HQET and read (Eq. (12) of Roca, Oset and Singh [22 



1-Vt°(s) G l=0 (s) 



(37) 



This equation manifestly is a rclativistic generalization of 
the Lippmann-Schwingcr equation. 

The factorized resolvent function is the standard one- 
loop integral 



4.2 Unitarized scattering amplitude 

Chiral perturbation theory amplitudes are by construc- 
tion a series expansion (albeit with logarithmic corrections 



Gi =0 (s) = i 



d 4 q 



1 



(2tt)4 (P - qf 



Ml 



it q z 



(38) 
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We employ dimensional regularization of the divergent in- Heavy-quark spin symmetry dictates that, whether D or 



tegral to read (Eq. (14) from ref. [22]) 



16tt 2 



a{n) 



M 2 
ln^% 
M 



Ml 



2s 



G l=0 {s) = 
Ml 



(39) 



-ln(s 



-\ [ln( S - (Ml - ml) + 2qV~s) 
V s 

+ ln(s + (Ml - ml) + 2q^fs) 
+- ln(s - (Ml - ml) - 2qV~s) 



(M 2 D -ml)- 2qV~s) - 2m] 



where the imaginary part of the logarithms above Dtt 
threshold reads 



3 Gz=o(s) 



8ny/s ' 



(40) 



with q the modulus of the pion's three-momentum in the 
CM frame. 

Introducing the conventional two-body phase space 




this imaginary part is 



3 G l=0 (s) 



PirD(S) 

16tt 



(43) 



With these ingredients it is straightforward to show 
that, by construction, the complex T a 's satisfy single- 
channel unitarity relations 



ZTt°(s) = -\Tt { '( 



1=0, a) * 



Pttd(s) 
16tt 2 



(44) 



(providing a convenient numerical check of our computer 
programmes). The amplitude can be parametrized in terms 
of the phase-shift 



T^(s) 



sm5 I0 (s)e iSl0 ^ 



Pwd(s) 



that arc then extracted via 



tan 5 io (s) = 



5 T w (s) 

nr I0 (s) 



(45) 



(46) 



with / = 1/2,3/2. (The tangent extraction should auto- 
matically resolve the phase-shift sign) . Finally, the isospin 
averaged amplitude for the LO-HQET decoupled single- 
channel problem becomes 



|T a | 2 = i(2|T^0|2 +4 | T 3/2,0|2 



(47) 



D* in any spin state, the scattering cross-section will be 
the same, and since an s-wave cannot flip the spin upon 
interaction, no further spin averaging is needed in leading 
order HQET. One can then use 



J2\M nc (s,t, x )\ 2 = \T a f 



(48) 



in Eq. 

Going beyond LO in HQET we need to distinguish 
between Dtt —> Dtt and D*ir —> D*ir scattering. To im- 
plement it, we assume that a charm quark propagates as 
a linear combination of both states 



\c)=a\D) 



■P-\T>* 



(49) 



The moduli of the complex numbers a and /3.; are deter- 
mined by thermal Bosc-Einstcin distribution factors, since 
the mass difference between D and D* slightly suppresses 
the latter. We then average over the relative (quasi-random) 
phases of a and (3 upon squaring to construct ^ |A^ wc (s, t, x)| 2 - 

For ease of comparison with other systems, we will also 
quote numerical results for the cross-sections given by 



and 



a(s) 7ri 



a(s) n i 



1 



16-7TS 
1 



D 



16irs 



D" 



(50) 



(51) 



although what is substituted in the Fokker-Planck inte- 
grals is the squared matrix element of M. . 



4.3 Value of the low-energy constants 

In the philosophy of low-energy effective theories, after 
all the symmetries have been used to constrain the La- 
grangian density, the remaining free constants have to 
be fit to experimental data. Eventually these constants 
should also be accessible to lattice QCD. 

To the order that we are working, the pion decay con- 
stant in the chiral limit F can be approximated by its 
physical value, = 92 McV, the difference being of one 
higher order in the chiral expansion. 

The renormalization scale for the NLO ChPT con- 
stants is to be understood as n = 770 MeV, and the 
scheme is such that the subtraction constant a(/i) = 1.85 
is fixed as in Oset, Roca and Singh [2"2~| . 

The authors of Ref. [S] quote a value of g = 1177 ± 
137 MeV for the heavy-light pseudoscalar-vector coupling 
constant <?, that can be obtained by reproducing the decay 
of D* + -mesons. We reproduce this elementary calculation 
with the Lagrangian density in Eq. (|18[) and obtain 



r = g 2 \p^/(l27TF 2 Ml»] 



(52) 



in agreement with [5], whose value and error band we 
adopt. 

In his recent paper [3], Laine quotes the value g T ~ 0.5 
for his effective Lagrangian. This Lagrangian is worked 
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out in detail in the textbook of Manohar and Wise [25] . 
where they quote an early value of g-^ = 0.42 from a lattice 
Montecarlo simulation by the UKQCD Collaboration. 

However both references employ a representation based 
on a heavy-hadron spinor field H a with dimension 3/2, 
whereas in our Lagrangian the D-field's dimension is 1. 
For these reason their g v has no dimension whereas our g 
has dimension = 1. Direct comparison in the Lagrangian 
is not transparent, but instead one can easily compare the 
D* tree-level decay width, and find the relation among the 
two couplings. The decay width employing the convention 
of Manohar and Wise reads 

r = gl\ P ^/{^) 

yielding g = \/2g n MD* ~ 1190 MeV, in agreement with 
the value from [5]. 

Turning now to the NLO constants, we have repeat- 
edly stated that hi = hi is a requirement of heavy quark 
symmetry tying the D and D* amplitudes at LO in Heavy 
Quark Effective Theory. We, in this article, set there- 
fore hi — hi = from the start. Likewise we have set 
ho — hi = hi = based purely on large-iV c count- 
ing. These constants well deserve being revisited in future 
work, but we are content here with accepting a l/N c sys- 
tematic error as customary in the current literature. 

Another useful constraint is offered by the mass differ- 
ences between the D-mcsons 3,5 , which fixes hi rj —0.45. 
Thus, the remaining free LECs to be estimated are 
and h§. We have at our disposal, in the Dtt channel cor- 
responding to the T a scattering amplitude, two pieces of 
known data (the Dq mass and width) to which we can tie 
/13 and /15. If a calculation including the subleading order 
in HQET is performed, then the D\ and Do parameters 
differ and the constants become overconstrained by known 
data. We will find in section [S] that reasonable values are 
(h 3 , h 5 ) = (7 ± 2, -0.5 ± 0.2 GeV" 2 ) with correlated er- 
rors, that is, an increased /13 needs to be used with a more 
negative h§. 

A word of caution seems convenient about the numer- 
ical value of /13 and h$. In [4] it has been proposed that 
h 5 = (h' b /m 2 D ) ~ 0{l/m 2 D ) since h' 5 is assumed there to be 
of order 1 . However we think this is unnaturally small and 
that /i,5 = 0((? 2 //1q CD ) should be expected. Our reason- 
ing is based on resonance saturation. Instead of unitarizing 
the amplitude and fitting the constants to the dynamically 
generated resonances, we could have introduced the reso- 
nances as additional fields and eliminated them from the 
low-energy theory [18] by employing 

-i 

' p 2 - m 2 Da (m D +m 1T ) 2 - m 2 Do 

near threshold. It is clear that the denominator is propor- 
tional to the off-shellness of the Dq resonance and not to 
its total mass. Of course, the analogous quantities coin- 
cide in the traditional case of 7T7t scattering since pions 
are so light as compared to the p for example, = 138 
MeV < m p = 770 MeV. Then p n7T ~ 0. However the mass 
of the ground-state charmed meson cannot be neglected 



in D-pion scattering and the low-energy constants do not 
vanish in the heavy-quark limit. Thus we would expect 
the denominator to be of order ^1q CD or at most m^m^, 
but not m 2 D . Some additional discussion about the hi con- 
stants can be found in a recent paper in the heavy quark 
limit [12]. 

5 Numerical results 

5.1 Cross section for Dtt elastic scattering 

We now present numeric computations of the unitarized 
and squared amplitudes in Eq. (f3"T)) . and of the cross- 
section. 

In the first place, and to compare with the work of 
Gamermann and Oset [5], we keep only the (s — u) term 
in the Dtt elastic amplitude V a . The square amplitudes 
with isospin 7=1/2 and 1 = 3/2 and I = are depicted 
in Fig. [31 The figure shows how the exotic I = 3/2 is non- 
resonant (this will also be the case for all the calculations 
presented below), which could have been guessed because 
no qq state exists with such isospin, so there is no intrinsic 
strength at low energies in exotic waves. The non-exotic 
1=1/2 channel presents a clear s-wave resonance, with 
approximate mass and width M ~ 2140 MeV and r ~ 
170 MeV. These values are somewhat too low if compared 
with the experimental M Do = 2360(40) MeV and r Da = 
270(50) MeV taken from the Review of Particle Physics. 

We do not deem this a problem since there is room 
for the NLO terms containing the hi constants to modify 
the computation. But how can then Gamermann and Oset 
obtain reasonable agreement with the experimental state, 
employing only the leading order amplitude in ChPT? 

We believe to have identified the reason in their sub- 
stituting one of the powers of /„• by /d, 

C q (s - u) C (s - u) 

This suppresses the strength of V a such that T a saturates 
unitarity at a higher center of mass-energy y/s, in better 
agreement with experimental data. Reducing V a by a fac- 
tor 2 displaces the maximum of the cross section to about 
2320 MeV with width about 250 MeV. 

The substitution of / w by Jd can be tracked to those 
authors employing SU{4) symmetry to construct the ef- 
fective Lagrangian, treating D mesons on equal footing 
with pions. However we believe this is a questionable pro- 
cedure since SU(4) is not even an approximate symmetry, 
and we have instead constructed the chiral coupling of 
pions to the heavy D-meson source. 

The LO interaction is therefore somewhat too strong 
in our case. We could weaken it by modifying the sub- 
traction constant of the loop function a(fx), but instead 
we proceed to the next order in chiral perturbation the- 
ory, which should be equivalent as shifts in a(p.) should be 
absorbed in the NLO hi constants. 

Next we add one by one the NLO constants hi, hs, 

Because the hi term does not increase with momen- 
tum, but is multiplied by a small m 2 constant, it does not 
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Fig. 3. Top: square amplitudes for Dir s-wave elastic scatter- 
ing employing only the (s — u) term of the interaction poten- 
tial V a (as in Gamermann and Oset). Bottom: isospin averaged 
cross section associated to those amplitudes. 



change the amplitudes appreciably. We include it but do 
not discuss it any further. 

We examine then the sensitivity to h% in Fig. [4] For 
small, positive values of /13 the Dq peak moves to larger 
masses, with /13 = 2 the shift is of order 50 MeV. The 
resonance also becomes broader. Then, for larger values 
/13 > 3, the mass starts falling again, and a cancellation 
with the s — u term sets in, forcing a zero of the amplitude 
at energies 2.5 GeV or above. 

If we now add the /15 term, we observe that its presence 
(if the sign is chosen negative as in Guo et al., for example 
/15 = —0.25 GeV -2 ) narrows the resonance shifting it to 
slightly lower masses. If positive, h$ forces a cancellation 
(as did a large h 3 ) giving a zero near threshold (for h 5 ~ 1 
GeV -2 , /13 ~ 2) or at 2.4 GeV and above (for the same 
h 3 but h 5 ~ 0.25 - 0.5 GeV -2 ). 

Therefore a strategy to improve agreement with the 
experimental Dq data is to combine a positive /13 with a 
negative /15 to increase the resonance mass without dis- 
torting the line-shape unacccptably. 




E CM (GeV) 

Fig. 4. Squared isospin 1/2 amplitude for D-k scattering for 
various values of /13, from top to bottom being 0,1,2,3,4. In this 
graph hs is kept fixed at zero. 
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Fig. 5. Squared isospin 1/2 amplitude for Dn scattering with 
(h3,h 5 ) = (7, -0.5 GeV -2 ) (central value). 



Our best computation is then shown in Fig. [SJ 
The maximum of the squared amplitude, employing 
(h 3 ,h 5 ) = (7,-0.5 GeV -2 ) as central value, gives a rea- 
sonable Md = 2300 MeV, just slightly below the exper- 
imental value, and a width just slightly above r — 350 
MeV. The two parameters are very correlated, so that 
varying one significantly requires varying the other simul- 
taneously to maintain reasonable agreement with the ex- 
perimental resonance. Shown in the figure are two more 
lines with the error band Ah^ = ±2 and Ah$ = ±0.2 
GeV -2 . It is this squared amplitude, leading order in Heavy 
Quark Effective Theory, that we adopt in our Fokker- 
Planck equation for the transport coefficients. 

Although the diffusion and drag coefficients require the 
\A4\ 2 square amplitude, it is convenient for the discussion 
to also plot the resulting cross section, which we do in 
Fig. El 
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Cross section for Dn elastic scattering with (/13, /15) 
GeV - ) (central value). 
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Fig. 8. Effect of including the Born terms associated with the 
D* . The bottom line (purple) is the cross-section associated to 
the Born term alone, as in the model of [TT]. The top line (red 
squares) is the resulting cross-section combining the Born term 
with the contact terms, without modifying the hi constants 
from Fig. HJ] and then unitarizing. 
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Fig. 7. Effect of adding the hi counterterms to the s — u basic 
Dn amplitude for isospin 3/2. 



The maximum of the cross-section is about 13.5 ± 1 
mbarn, and for the entire range of center of mass energies 
y/s £ (2 — 3) GeV we find a > 5 mbarn. In fact, for a large 
window between 2.1 and 2.5 GeV we have a > 10 mbarn, 
which is slightly larger but in reasonable agreement with 
the guess by the authors of [10], that assume 7— 10 mbarn, 
or by Svetitsky and Uziel [T] of 9 mbarn. 

For the sake of completeness, we separately quote the 
effect of adding the hi constants on the non-resonant isospin 
3/2 Dtt elastic amplitude. The corresponding plot is num- 
ber 

As can be seen in the figure, the effect is moderate at 
all energies. 

Next, we proceed to the next-to-leading order in Heavy 
Quark Effective Theory. We only consider for now the 
Born s and i-channel exchange terms due to D* exchange 
between the Dtt pair. 

The effect of adding these terms is akin to making /15 
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Fig. 9. Effect of including the Born terms associated with the 
D* , but leaving the hi coefficients free. The red, solid line is 
the central value with /13 = 8, /15 = 0.35 GeV~ 2 . The black, 
dashed line coincides with the cross section in Fig. [15] without 
the Born terms. 



more negative, that is, a narrowing of the Dq resonance, 
as shown in Fig. [5] 

However, a renormalization of the hi constants effec- 
tively brings back the pole position in better agreement 
with experimental data. We now refer the reader to Fig. [9] 
Shown in the figure are lines with (/13, /15) = (7.5±2.5, 0.4± 
0.3 GeV -2 ), together with the result of Fig. [B]without in- 
cluding the Born terms. As can be seen, the effect of the 
D* exchanges can be largely absorbed in the hi countert- 
erms (for fixed m c mass of course, since they scale differ- 
ently) so we will ignore the Born terms in this computa- 
tion. However a certain uncertainty should be understood, 
of order 30% in the cross section , that could be larger than 
our estimate in the region of the Dq resonance. 

It is also worth commenting that the addition of the 
Born terms causes a dip in the high-mass Dtt spectrum 
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Fig. 10. Elastic cross-section for D*n scattering computed 
replacing mo by m D in Fig. [6] The resonance should now be 
interpreted as the broad Di(2427). 



that can be brought down by minimum changes in the hi 
parameters. Since we do not think that such an interfer- 
ence dip between D* Born exchange and the Dq tail has 
been reported in experimental data, we keep the contact 
parameters in a band such as not seeing this dip in the 
momentum range of relevance. 

Finally, we return to the computation in Fig. [6l but 
substitute mn by m* D (an NLO effect in HQET) as only 
modification to obtain Vd instead of V a ■ We interpret the 
resulting cross-section as that corresponding to D*ir scat- 
tering, and plot the result in Fig. [TO] The cross-section 
including both 1/2 and 3/2 isospin channels is clearly res- 
onant, with the D\ well visible. As was the case for the 
Do, the mass is slightly below the data. The cross-section 
peak is about 15 mbarn. 

Thus we have performed an exhaustive study of the 
LO-HQET interaction and now proceed to compute trans- 
port coefficients equipped with the interaction leading to 
Figs. [U and [TO] 



5.2 Diffusion and drag coefficients 

We now proceed to computing, with the square ampli- 
tude so numerically computed, the F, Jq and i~i transport 
coefficients. The three pannels of Fig. [TT] shows them as 
function of squared momentum p 2 for fixed temperature 
T = 150 MeV. One should not trust these results above 
charm momenta of order p = 1.5 GeV, but we spell them 
out for completeness. 

In the top panel of this figure we show the drag coef- 
ficient F(p) in fm _1 , which exhibits a momentum depen- 
dence of about 10 % within the range of p G (0, 2.5) GeV. 
From this coefficient one can extract the relaxation length 
for a charm quark propagating in the pion medium that 
turns out to be around 40 fm at p =1 GeV. 

Quite strikingly, one can see in the figure that Jq has a 
very mild momentum dependence, its value can very well 
be approximated by a constant for the entire momentum 
range. J\ is seen to grow with momentum, increasing the 
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Fig. 11. We show all three coefficients in the Fokker-Planck 
equation as a function of charm-quark momentum, at a refer- 
ence temperature of 150 MeV in the pion gas. The low-energy 
constants in the Dir — ► Dir amplitude are fixed to hi = —0.45, 
g — 1177 MeV, and /13 and /15 fit to describe the mass and 
width of the Do resonance. Top: F including and not including 
the possible propagation of the c quark as a D* meson. Middle: 
7~b and A including D-like propagation alone. Bottom: Fo and 
Fi including also propagation as a D* meson. 
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Fig. 12. Momentum-space drag coefficient as function of tem- 
perature for a stopped charm quark in the hadron gas. We ob- 
tain the coefficient by employing the Einstein relation in taking 
the limit of p — > 0. 



difference Pi — To, and thus favoring diffusion at higher 
typical momenta. 

In Figs. [T^] and Q2] we show the dependence with the 
temperature of the drag coefficient at fixed momentum. 
Since the direct computation of F(p 2 — > 0) is rather unsta- 
ble, the plot in Fig. [T^] is computed from F by employing 
the Einstein relation, Eq. (|103[) in the Appendix. 

The drag coefficient is seen to increase by a factor of 
about 4 in the range from 100 to 150 MeV, so that most of 
the drag in a heavy ion collision is expected in the hotter 
stages, with the charm quarks freezing out progressively 
until they freely stream outwards till they decay. 

We compare with other authors, choosing a reference 
temperature of 100 MeV where all existing works make a 
statement, and show the drag coefficient for each recent 
work in Table [3l It can be seen that the phenomenological 



Table 3. Value of the drag coefficient at p 
MeV. 



and T = 100 



Authors 


F(fm^ i ) 


Laine 


0.05 x 10~ 3 


He, Fries, Rapp 
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Ghosh et al. 
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0.04 



£ 0.03 



0.02 
0.01 



CM 
> 

CD 

o 



> 


o 



0.016 
0.014 
0.012 
0.01 
0.008 
0.006 
0.004 
0.002 




0.024 
0.022 

0.02 
0.018 
0.016 
0.014 
0.012 

0.01 
0.008 
0.006 
0.004 
0.002 




— p=1 GeV 

p=0.6 GeV 

p=0.3 GeV 
p=0.1 GeV 



100 110 120 130 140 150 160 170 180 

T (MeV) 



p=1 GeV 
p=0.6 GeV 
p=0.3 GeV 
p=0.1 GeV 




100 110 120 130 140 150 160 170 180 

T (MeV) 



— p=1 GeV 

p=0.6 GeV 

p=0.3 GeV 
■ - ■ p=0.1 GeV 




100 110 120 130 140 150 160 170 

T (MeV) 



180 



Fig. 13. Momentum-space drag and diffusion coefficients as 
function of temperature for a slow charm quark with momen- 
tum p = 1 GeV, 0.6 GeV, 0.3 GeV and 0.1 GeV. Note that 
the intensity of the drag force is roughly proportional to the 
temperature. 



model of He, Fries and Rapp is of the same order of magni- 
tude of our result, with Ghosh et al. quoting an extremely 
large value in their Fig. 2, and Laine a smaller value by 
one order of magnitude. We believe that we have a larger 
control of the charm-pion scattering amplitudes at mod- 
erate temperatures, but the reader would be cautious to 
employ a factor 2 as error band to our result. 

The spatial diffusion coefficient is then plotted in Fig.[T4l 
as function of temperature. 



At low temperatures it correctly takes the non-relativistic 
limit 

3T 3/2 

D x = —= (53) 

with m the particle mass, a the cross section, and P the 
pion gas pressure, that is temperature dependent. We also 
note that, during the lifetime of the pion gas after the 
cross-over from the quark-gluon plasma phase, the inter- 
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Fig. 14. Spatial diffusion coefficient as a function of temper- 
ature. 



100 
90 
80 
70 
60 
50 
40 
30 



20 
10 



H =-100 MeV 
H*=-50 MeV 
H*=0 MeV 
H*=50 MeV 
n"=100 MeV 
n"=138 MeV 



100 110 120 130 140 150 160 170 180 

T (MeV) 

Fig. 15. Same as in Fig. 1 141 but as a function of the chemical 
potential 

actions between pions are almost entirely elastic, so that 
pion number is effectively conserved and one should intro- 
duce a pion chemical potential, not included in the very 
recent works by other groups. Introducing this approxi- 
mate pion chemical potential [i, 



P oc mi /2 T 5 ^ 



(54) 



makes the product TD X diverge at low temperature and 
vanishing chemical potential (which just means that gas 
particles are too cold and slow to stop the charm quark 
from diffusively moving inside the pion gas). However, at 
chemical equilibrium with (i — > (that is not expected 
in the hadron phase of a heavy-ion collision, but is rele- 
vant to make contact with the non-relativistic limit), the 
exponential becomes unity and TD X becomes a constant 
at low temperature. We further show the effect of this pion 
chemical potential in Fig. [15j 

We find the effect sizeable. At a reference temperature 
of 120 MeV, the ratio between D x at /!„■ = and \i v = 138 
MeV is a factor of about 5. 

To assist in the physical interpretation of these results, 
we have plotted in Figs. [Tfjl and [T71 the loss of energy and 
momentum per unit length discussed in Appendix [Cl and 
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Fig. 16. Loss of energy of a charmed meson as function of 
the energy in a pion gas at a fixed temperature of 150 MeV, 
assuming it can travel as a D or a D* meson during the few 
fermi of the gas's lifetime. 
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Fig. 17. Loss of momentum per unit length as function of mo- 
mentum of a charmed meson in a pion gas, same as in Fig. [T75J 



derived from our results for the drag coefficient F, for 
various momenta p 2 . 

From Fig. 1171 one can estimate that a reference charm 
quark in a D or D* meson with momentum 1 GeV mea- 
sured in the rest frame of the pion fluid surrounding it, 
will deposit about 50 MeV per Fermi travelled in the fluid. 
Thus, if the pion gas is in existence for, say, 4 fm, the D 
meson measured in the final state with a momentum of 800 
MeV will have been emitted from the quark-gluon plasma 
phase with a GeV. This result is similar to the 20% effect 
recently quoted by He, Fries and Rapp [TU] and means 
that, while the D and D* mesons can be used as probes 
of the quark-gluon plasma, their distributions should be 
shifted up in momentum (or alternatively both the quark- 
gluon plasma and hadron phases have to be treated in 
hydrodynamic simulations). 

The authors of reference [12] proposed to divide the 
temperature times the spatial diffusion coefficient by the 
shear viscosity over entropy density ratio ry/s, produc- 
ing a dimensionless quantity that should give an idea of 
how strongly coupled is the quark-gluon plasma, and they 
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Fig. 18. A dimensionless ratio with the viscosity over entropy 
density, proposed in Ref. [13]. The top dashed line corresponds 
to a "weakly coupled" quark-gluon plasma, the bottom line 
to a "strongly coupled" quark-gluon plasma. The solid line at 
the bottom, for charm propagating in our pion gas, is more 
suggestive of the second than of the first. 



quote two estimates based on AdS-CFT that we plot in 
Fig. [THJ In the figure we also plot our computation based 
on charm quarks travelling through the pion gas, together 
with our computation of viscosity over entropy density in 
the pion gas presented in [28 a . It seems that, according 
to this criterion, the charm quark is somewhat strongly 
coupled to the pion gas, although it is not clear what the 
precise value of these AdS-CFT based estimates is. 



6 Experimental discussion 

While we are not directly computing the experimentally 
observed quantities in this article, it is worth looking ahead 
onto what impact our results have for the Heavy-Ion col- 
lision programme at RHIC and at the LHC. 

A commonly quoted observable is the nuclear suppres- 
sion factor Raa obtained by dividing the number of elec- 
trons from heavy meson decays in a nucleus-nucleus colli- 
sion by the number in proton-proton collisions times the 
number of constituent nucleons, 



Raj 



N aa 



A 2 x N, 



VP 



At high pr of order 5-7 GcV this ratio reaches 0.3, show- 
ing substantial effects due to the medium. At low trans- 
verse momentum up to to 2 GcV this suppression factor is 
close to 1 (small effect). At the lowest pr of few hundred 
MeV the ratio is even larger than unity, there being an 
enhancement of the number of heavy mesons in ion-ion 
collisions [23] over proton-proton. This pile-up of heavy 
mesons at low momentum can be interpreted as they being 
slowed down by the medium. The F and r± coefficients are 
relevant for this process, with J\ broadening the px dis- 
tribution and F equating the velocity of the heavy quarks 
to the velocity of the fluid medium that they are crossing, 
as can be seen in Eq. ([5"5]) below. 



Another important observable is the elliptic flow [25] 
V2 defined by the distribution of particles with the az- 
imuthal angle <f> taken around the collision axis, with the 
collision plane at = 0, 



1 dN 



1 



N dej) 2n 



(1 + 2U\ cos <j) + 2z*2 cos 2(j) + 



Substantial elliptic flow for heavy flavored mesons has 
been measured at RHIC [23], meaning that the heavy 
quarks are partly equilibrating with the medium. This el- 
liptic flow can potentially provide sensitivity to the com- 
bination of diffusion coefficients J\ —Tq. In subsection !5.2l 
we showed that this difference grows with quark momen- 
tum, which may help explain why the elliptic flow grows 
in the pr 0-2 GeV range. 

As for the actual spectrum of D-mesons deduced by 
STAR [25] , although the errors are very large, they quote 
an average velocity of j3 = 0.35 — 0.47, sufficiently smaller 
than 1 to make HQET a reasonable starting point, espe- 
cially taking into account that part of this velocity is due 
to the local fluid rest frame being in motion in the labora- 
tory frame, with the actual charm velocity respect to that 
Eulerian frame being even smaller. The STAR collabora- 
tion also quotes a rough temperature of 120 MeV as fitting 
their spectrum, but given that they have only three points 
in the plot and the large error bars, a larger temperature 
(or a poor thermal fit) would not be surprising at all. We 
should wait for future data to clarify this point. 

These observables have not yet been provided by di- 
rect reconstruction of the D or B mesons (maybe AL- 
ICE can provide a measurement with its Internal Track- 
ing System assisting the secondary vertex reconstruction), 
but indirectly with measurements of the (presumed) sec- 
ondary electrons from heavy flavor decays. Copious D- 
meson counts have already been informally reported in 
the Kir and Ktttt channel and we look forward to the 
publication of this data. 

However we feel that the measurements are very en- 
couraging and that we should expect these charm drag 
and diffusion coefficients to become accessible. Then it 
will be necessary to disentangle the conventional diffusion 
in the hadron phase from the more exotic quark and gluon 
phases, and our results will be useful here insofar as they 
greatly reduce the uncertainty in the hadron, low temper- 
ature phase. 

In Fig. [201 of the appendix we show the solution to 
a one-dimensional version of the Fokkcr-Planck equation 
that comes handy for this discussion (assuming mid-rapidity 
and no azimuthal flow, the charm diffusion, although more 
complicated, is reminiscent of that one-dimensional case) . 

From that figure one can see that if the initial distribu- 
tion of charm quarks would peak at some po € (1, 2) GeV, 
for every femtometer spent in the pion gas, the charm 
quark distribution would peak 50 MeV lower (friction), 
and the distribution would be about 100 MeV broader 
(diffusion) . 

In fact, the ALICE collaboration has already published 
an analysis for proton-proton collisions to serve as bench- 
mark |23 for what is to come in Pb-Pb. As usual, the mid- 
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Fig. 19. We compare our computation of the diffusion co- 
efficient with other estimates. The possibility of a minimum 
of the charm relaxation time at the phase transition seems to 
be well possible. The leftmost dashed line is pure perturbation 
theory [9], to which our result seems to tend asymptotically 
at low temperature. In addition, we have plotted two curves 
below and above the crossover [101115] 



rapidity p t D-mcson spectrum has an exponential shape 
~dN 



dp t 



oc e A 



(55) 



y<0.5 



with a scale A ~ 1.5 GeV and it will be interesting to 
convolute our Fokker-Planck kernel with this input exper- 
imental spectrum. 



7 Summary and conclusions 

If we take 150 MeV as the highest temperature at which 
our approach is reliable (as we do not include strange 
mesons), and we examine a charm quark travelling as a 
meson with momentum 1 GeV, the relaxation length read 
off from the top plot in Fig. [TT] is 



X c (p = 1 GcV,T = 150 MeV) 



1 



0.025 fm 



40 fm 



happens around the phase crossover, where the interac- 
tions also have longest range and intensity. Thus charm 
quarks can be potentially used as a probe of the phase 
transition if theoretical uncertainties on the hadron gas 
side can be reduced. We believe that we have produced 
a very reliable estimate of the hadron coefficients in the 
temperature region T < 150 MeVQ 

Laine finds a formula for the momentum space diffu- 
sion coefficient 



This is much larger than the expected lifetime of the hadron 
gas A ~ 5 — 10 fm. He, Fries and Rapp have reported re- 
laxation times as short as 25-40 fm in the 150-180 MeV 
region, in good agreement with our estimate. It should 
be taken into account that the hadron gas cools down to 
lower temperatures in the last stages, with longer relax- 
ation times. Thus, charm quarks will not completely relax 
during the lifetime of the pion gas and will indeed carry in- 
formation from the crossover from the quark-gluon plasma 
phase. 

In Fig. [T5] we show a comparison of our computation 
with the best existing ones and with the estimate of [T5] 
that employs a mixed plasma plus resonance approach 
above the crossover. All existing information points out 
to that the minimum relaxation time of the charm quark 



r oc 



T_ 

f 4 



provided that m^/n <C T <C fn, which is a very restrictive 
range of temperatures around 60 MeV. We have shown 
that this growth with temperature is way too fast and that 
properly unitarizing the interaction tames this high power 
of the temperature. We also qualify the statement that, in 
this range, the coefficients are dependent only on the pion 
mass and decay constant; this should be understood as 
valid in the infinite quark mass limit, while in the charm 
sector we find that the m* D — mass difference brings 
about non-negligible corrections. 

Svetitsky and Uziel (Fig. 1 in [1]) found that a c-quark 
with initial transverse momentum 2 GeV would have come 
down to 1 GeV by the time of freeze out. What our re- 
sults show is that all this decrease needs to be assigned 
to the quark and gluon plasma phase and, especially, to 
the phase transition, but that the loss of momentum in 
the pion gas is a moderate-sized correction. For example, 
a 1 GeV charm quark entering the pion gas at 150 MeV 
and travelling four femtometers through it, will have lost 
about 200 MeV at freeze-out. 

We have found that the F drag coefficient and the 
To diffusion coefficient depend only mildly on the charm- 
quark momentum, implying that the nuclear suppression 
factor for charm in the pion gas can be reasonably mod- 
elled. On the contrary, we find that the A diffusion coeffi- 
cient strongly depends on momentum, so that anisotropic 
observables such as the elliptic flow will have a more in- 
volved dependence with momentum. 

Moreover we have shown that the thermal relaxation 
time at 150 MeV is about 40 fm, implying that the charm 
quarks do carry memory of the phase transition upon ex- 
iting the hadron gas. Our results also suggest that the 
spatial diffusion coefficient is likely to have a minimum at 
the crossover to the quark and gluon plasma. 



1 Our numerical data for scattering amplitudes or transport 
coefficients is at the disposal of interested colleagues who want 
to pursue kinetic or hydrodynamic simulations by contacting 
fllanes@fis.ucm.es or j .torres@fis.ucm.es. 
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A Static Fokker-Planck equation and 
one-dimensional solution 

Suppose that the three scalar coefficients F, Jo, A in 
Eq. (|14p do not depend on p (limit of momentum inde- 
pendence or static limit). Then one can speak of constant 
F{p 2 ) = F (we will now show that this is simply a friction 
coefficient), and r (p 2 ) = A(p 2 ) = r (diffusion coeffi- 
cient). 

The Fokker-Planck equation reads 



^ = F\7 p - (p/ c ) + r\/ 2 p f c 



(56) 



that can be compared with the standard diffusion equation 
for the concentration C of a solute 



^ = - A A7-(CF)+#V 2 C, 



where (i is the mobility, D is the diffusion coefficient and 
F is the external force. Einstein's relation D = fiT relates 
the diffusion coefficient and the mobility. 

For clarity let us concentrate on one dimension. In this 
simpler case, the equation 



— = F—( f ) [ 

dt dp c dp 2 



(58) 



is known as Raylcigh's equation and describes the momen- 
tum distribution equation for a Brownian particle. 
With the initial condition 



f c (p,t = 0) = S(p- Pa ) 
the analytic solution reads 



(59) 



fc(p,t) = 



F 
2^T 



l-e 



-2Ft\ 



-1/2 



exp 



F (p-poe 



-Ft\2 



2F 1 



,-2Ft 



(60) 
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Fig. 20. One-dimensional evolution of a drop of charm f(t,p) 
from Eq. (|60|l with well-defined momentum po = 1, 1.5 and 2 
GeV (top to bottom) , at a reference temperature of 150 MeV. 
Time evolution drags the momentum towards zero from the 
initial condition, Dirac's delta function in Eq. (|59[) . and the 
shape broadens to adopt the Boltzmann equilibrium function. 



We can identify the large-time behavior of this function 
with the Maxwell-Boltzmann equilibrium function 



lim f(t,p) = /mb(p), 



(61) 



This solution can be easily plotted, and we do so in Figure provided that an analogous to Einstein's relation holds 
[2D] for momentum po = 1, 1.5 and 2 GeV (top to bottom), 

at a reference temperature of 150 MeV. r = FmT, (62) 
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and the two coefficients F, r are not independent but 
rather related by this fluctuation-dissipation relation (as 
shown below in Eq. (|103[) where we have derived this re- 
lation from the momentum-dependent fluctuation dissi- 
pation theorem independently of the number of spatial 
dimensions). Morever, in Appendix [Bl we show that the 
momentum diffusion coefficient, r, is related to the spa- 
tial diffusion coefficient, D x , as: 



D, 



r 



2 p2 



T 2 

T' 



(63) 



B Classical Langevin equation for charm 
diffusion 

The purpose of this appendix is to show the relation be- 
tween the diffusion coefficient in space D x (that appears 
in Fick's diffusion Law) and the momentum diffusion co- 
efficient D, that we have estimated through the Fokker- 
Planck equation. This discussion is well-known from clas- 
sical statistical physics, but it is enlightening to review 
it and makes the article self-contained. Since we use this 
material mostly to give the various quantities a physical 
interpretation, we believe that it is sufficient to limit our- 
selves to a purely classical discussion (as appropriate for 
a dilute gas). 

We begin by rederiving the Fokker-Planck equation 
from the Langevin equations. Several manipulations of 
Dirac delta distributions are easier to follow discretizing 
the time variable, to avoid resource to somewhat advanced 
functional analysis, and we will thereafter take again the 
continuum limit St — > 0. Then the classical solution to the 
Langevin will allow us to identify the space-diffusion term 
and relate it to the Fokker Planck coefficient of diffusion 
in momentum space. 

The charm quark (Brownian particle) moves in the 
pion gas and it is diffused because of collisions with these 
mesons. The position and momentum of the charm quark 
can be regarded as stochastic variables depending on time. 
The classical, non-relativistic stochastic differential equa- 
tions that govern their motion are: 



dx 1 p l 
dt mo 
dp 1 
~dt 



F l (p)+e(t) 



(64) 
(65) 



where the index i = 1,2,3 labels the space component of 
x and p. This equation is called the Langevin equation. 
The F z (p) is a deterministic drag force which depends on 
momentum through the collision processes and is a 
stochastic term called white noise. It verifies 

(£*(*)> =0, (66) 

(C(t)e(t')} = r^(p)S(t-t'). (67) 

In an isotropic gas one has -T y (p) = r(p)5^ . 
We now discretize the time variable 

t„ = nSt; x„ = x(i„); p„ = p(t„); n = 0, 1, 2, 



and choose a mid-point discretization for F [5D] 

^(p) = F'[ P " + 2 P " +1 ]- 
The discrctized Langevin equation reads then 

Pn 



(69) 



x„ + i = x„ + - — St , 



(70) 

Pn+1 = Pn - F„ 5t + L n St , (71) 
with a time average over the random noise 

Li = jJ* n+l dte{t). (72) 

From jBBJ and ([57]), L\ verifies: 

(L^ = . (73) 

(LIA) = §<^'<W • (74) 

(With somewhat more work one can show that the vari- 
able L l n ~ 0(St~ 1 / 2 ) ). 

The average () is taken with respect to the probability 
associated with the stochastic process. Since the stochas- 
tic variables are positions and momenta, this probability 
is nothing but the one-particle classical distribution func- 
tion, f(t, x, p). Averages arc then computed by means of 

(T(t)) x ,p = J dxdpT(t, ) , (75) 

where T(t, x n ,p n ) is any function of the stochastic vari- 
ables and time. 

In the Fokker-Planck equation we look for the time 
evolution of the distribution function itself, so we need to 
calculate the probability that a particle at time f n +i is at 
x,p 

/(t„ +1 ,x,p) = (<5( 3 >(x„ +1 -x)^( Pn+1 -p)) , (76) 

from the distribution function at a prior time. 

We introduce the discretized Langevin equation inside 
the deltas in (fTB")) : 



S(x n+ i - x) = S(Xn - x + —^St) 

m D 



(77) 



5{p n+x -p) = 5{p n P + [F„ + L n ] St) . (78) 
Expanding the deltas up to 0(St), 

-x*) = 

5(< - x l ) + £ ^jS(xi - x l ) ^St , (79) 
2 dx J n rn-D 

s( P l n+ i-p l ) = 

■ Op n 



(68) 



4ee 



-r 1 k dpkdpk 



Sipi-p^LiLl {5tf , (80) 
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and introducing these expansions inside equation (|76[) , we and replace the definition in Eq. (|76[) 
see that 



/(i n+ i,x,p) = (<5( 3 >(x„ - x)5( 3 )(p„ - p)> 



(5& (x„ - z)^ 3 ) (p„ - q)) = /(*„, z, q) . (83) 
One obtains 



+<£ A 5(3) ( x « - *) 5(3) (p« - p)>— 

j dx£ TO n (84) 

(^ 3 )(x„-x)V — <5 (3) (Pn-p) [4-^(p n )])5i /(W,X,p)= / dzdq^z-x^q-p) /(i n ,z,q) + 

5 9ri /• rl St 

j dzdq j( 3 )( q -p)^^L<5( 3 )( z -xV /(i n ,z,q) 



-i(<5< 3 )(x„ - x) lJ^ 6i3) (v- v)LiL k n )(5t) 

j k 



(gl) J dzrfq<5( 3 )(z-x)^A (5 (3) (q _ p)F , (q) /(tnjZiq)at + 

i 

dzdq <5 (3) (z - x) ^ — t ^-—<5 (3) (q - p) — - — f(t n ,z,q)St 

niuiuuuue une luiiuwuig iuenwu,y ij * * 

.( 3 \. \,r(3)/ \ / fi0 \ where the average operation has been factorized because 

{ n x j IPn PJ —I j om y depend on U n , with n' < n. 

J dzdqS^ (x„ - z)<^ 3 > (z - x)<J< 3 > (p„ - q)<5 3 (q - p) Now inte g rate b y P arts and finall y> over z and * 



/(x, p, t n+1 ) = /(i„,x, p) - -E- • ^/(*n,x, p)<ft + £ A F i( p ) /(ini X) p ) 5t + i ^ _JL_ r «(p)/(t„, x, p)*t . 



Now we can return to the continuum limit <5t — > 0: 



TO£> 9x" 



2 ^-r* dp 1 dpi 



r 



(85) 



Taking the average in space 

0/c(*,P) 



|"[^(p)/ c (i,p)]+i^-^(p)/ c( t, p) . 

(86) 



that coincides with the Fokker-Planck equation in Eq. (|T1?)) 
Now we see that the diffusion coefficients To, Pi, stem from 
the random force in the Langevin equations, and the drag 
coefficient from the deterministic friction force there. 

In the static limit p — > 0, we can solve the Langevin ^ s then solved by 
(or, in this limit, also Uhlenbeck-Orstein) equation 

dp 



fluid's rest frame. 

(p(t)) = Pa e- Ft . 
The second of Hamilton's equations 



dx 
~dt 



P 

m D ' 



dt 



(87) 



x(i) = x + / dr 
Jo 



P(r) 



whose solution is 



Or, on average, 



p(t) = p e 



-Ft 



-Ft 



dre FT £,{T) 



(x(t)> = x + 



Po 

Fmj 



■(l-e-«) 



(89) 



(90) 



(91) 



(92) 



Taking the average one can see that due to the drag force, 



To make the connection with the spatial diffussion co- 



thc friction term makes the particle eventually stop in the efficient we can show the mean quadratic displacement of 
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the Brownian particle (r = \/ x 2 + y 2 + z 2 ) 

((r(t) - r ) 2 ) = ((x(t) - .t ) 2 + (y(t) - y ) 2 + (z(t) - z ) 2 ) , 

(93) 

that, from Fick's diffusion law, is simply 

((r(t) - r ) 2 ) = 6DJ . (94) 
From the averaged solution to Langevin's equation (j9"T]) , 

<(^)-*o) 2 H-Jr 

m D JO JO 

With the help of and (|57|) and carefully perform- 
ing the integral [20] one obtains the leading term of this 
expression when t F^ 1 as 



drdr' (p x {t) P x {t 1 )) . (95) 



so that 



r 



m 2 D F 2 



2rt 

m 2 D F 2 



T 

T 1 



(96) 



(97) 



where finally we have used Einstein's relation. Thus, the 
calculation of the momentum space diffusion coefficient 
automatically entails an estimate for the space diffusion 
coefficient. 



C Fluctuation-Dissipation Relations 

Not all three coefficients F(p 2 ), r (p 2 ) and A(p 2 ) ap- 
pearing in the Fokkcr-Planck equation are independent, 
but rather related by a fluctuation-dissipation relation. 
This just means that thermal equilibrium requires the 
damping force F to match the fluctuations of the charm 
quark momentum distribution so as to maintain energy 
equipartition, with per degree of freedom. Since we 
consider the p-dependence of the three coefficients, the 
fluctuation-dissipation relation will be momentum depen- 
dent, although we also expose the p — > limit. More de- 
tails on deriving such relations can be found in standard 
textbooks [21j . 

A transparent procedure is to match the asymptotic 
solution of the Fokker-Planck equation to the thermal 
equilibrium distribution function, thus guaranteeing en- 
ergy equipartition. 

First of all, the Fokker-Planck equation can be written 
as an equation of continuity 



where 



df c (t,p) = _8_ 
dt dpi 1 







(98) 



n t = -Fi(p 2 )f c (t,p) - — [j^-(p 2 )/ c (*,P)] (99) 

is the particle flux density in momentum space. At sta- 
tistical equilibrium, this flux is zero, and the equilibrium 
distribution function is the Bosc-Einstcin function, 



fc 



1 

-p 2 /2AIT _ i 



(100) 
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Fig. 21. Momentum-space transport coefficients as function 
of temperature. Dotted: T (p 2 0). Dashed: A(p 2 -> 0). The 
very good agreement in our computer programme, as appro- 
priate in this limit, makes the curves barely distinguishable. 
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Fig. 22. Momentum-space transport coefficients as function of 
momentum at fixed temperature 150 MeV. The two coefficients 
converge at low p. 



Employing again the approximation 1 + / c « 1 valid 
for small charm-quark number, one can obtain 



dr^jp 2 ) 



(101) 



This momentum-dependent fluctuation-dissipation relation 
can be recast for the functions F(p 2 ), -To(p 2 ) and A(p 2 ) 
as: 



2 , 2 r 2 21 A(P 2 ) 

{p } + p dp + p 2 ' t l{p } _ o(p ^ " 



MT 



(102) 

For low-momentum charm quarks, /^(p ),Po(p ) — ^ 
r, F(p 2 ) — > F. The equality of the two r coefficients 
in the limit of zero momentum is numerically checked in 
Figs. EU and [22] 

We then recover the well known Einstein relationship 



F 



r 

MT 



(103) 
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which is the same result derived above in Appendix 1X1 for 
the one-dimensional solution of the Rayleigh equation. 

Thus, in the static limit two coefficients take the same 
value and the third is obtained from them by Eq. (|i03ll . 
and we are left with only one independent diffusion coef- 
ficient. 

The Langevin equation also allows us to directly obtain 
the classical interpretation of F as a loss of energy per unit 
length. Ignoring the fluctuating force, 



scattering amplitude to correctly position the Do and D\ 
resonances. For example, the outgoing-particle energies 
are E° 



k)2and^ +fe 



y/ml + (q + kf 



- P -k = V m D + 
respectively. 

The resulting transport coefficients obtained after inte- 
grating Eq. (|I07[) depend only on the modulus of p. How- 
ever we will introduce a (trivial) dfl p angular integration 
in the p-coordinates 



d-ymv 
dt 



(104) 



dk dq 



dkdq^ 

47T 



(110) 



can be multiplied by v to yield after some reshuffling the 
obvious expression for the power 



dt 



F v 



(105) 



and remembering the definition F = Fp in Eq. (Ti"4"j) . the 
loss of energy per unit length is simply F\p\, as in the 
non-rclativistic theory. 

The loss of momentum per unit length can then be ex- 
pressed as 

d ^ = ^r = -FE (106) 
dx vdt 

in terms of the energy and momentum of the charmed 
particle. 



D Kinematics 

In the evaluation of the drag and diffusion coefficients in 
equation (p~7|) we need to calculate integrals of the generic 
type 



in order to increase our freedom in the choice of axes. 

We also find convenient to change the integration vari- 
ables from the incoming pion and transfered momenta, q 
and k respectively, to the total momentum P = p + q and 
the outgoing charm momentum P3 = p — k. The Jacobian 
determinant associated to these translations is unity. 

Now, without loss of generality, we choose the total 
momentum P vector along the OZ axis, and the incom- 
ing charm momentum p lying on the OZX plane. Auto- 
matically q is in this plane as well. Finally, the outgoing 
charm momentum P3 has in general all three Cartesian 
projections, 

P = (0,0,P) (111) 
p=(py/l-xl,0,px p ) (112) 



q = -p-p=(-p^l-x%0,P-px p ) (113) 
^ = (y 1 - x l cos</>3, sin 03, x 3 ) . (114) 



97 



dkJ^Uq) [l + / w (q+k)] 



1 



1 



1 



(27r)9 j '< \ ~±/ v , . -„ ^ 2E c_ k 

(2n)*5(E c p + E* - E c p _ k - ^ +fc )E ^(M.xJlW.P*) 

(107) 

where the various cases differ in the choice of Q kinematic 
function 



kjp 1 
k 2 _ (fc.p 1 ) 2 



1 (fc»P*) 2 

2 p 2 



for F(p 2 ) 
forToCp 2 ) 
forA^ 2 ) 

The collision momenta are labelled as 

c(p) + 7r(q) — !• c(p - k) + 7r(q + k) 



(108) 



(109) 



Here x v is the cosine of the polar angle of p, that is, of 
the relative angle between p and P; x 3 and <fi 3 the cosine 
of the polar angle and azimuthal angle associated with 
p 3 . The transfered k = p p 3 and outgoing pion p\ = 
q + k = P — P3 momenta are then dependent variables. 

The angular integrals associated with P are then triv- 
ial (the scattering matrix is rotation-invariance) and yield 
47T (they have de facto being exchanged for df2 p that is 
now non-trivial). The system has one more rotational in- 
variance, as holding the P axis fixed, one can rigidly rotate 
all other vectors around it (so our choice of OX axis 
docs not imply any loss of generality) . This trivializes the 
<p p integration. 

With such choice of axes the integration measure can 
be explicitly written down as 



(115) 



so that = (Ep,p) is the 4-momentum of the incoming 
c quark, = (E* , q) the 4-momentum of the incoming 
pion, and k the transfered momentum from the c quark 
to the pion. The c quark can be in a D or in a D* meson 
states, degenerate in Leading Order Heavy Quark Effec- 
tive Theory. We generically use an average tod for the 
transport code, although we distinguish the masses in the 



dk dq — v - = / dP dp 3 — ^ 

47T J 47T 

4:TrP 2 dP p 2 dp3d(f>3dx3 — 2irdx p 

47T 

= 2.JdPP 2 d P3 pld^dx 3 dx p 



Energy conservation imposes an additional restriction, 
which is very non-linear in terms of the momentum vari- 
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K = \l m2 D+p 2 



ablcs 

(116) 

Eq = \]ml + P 2 +p 2 - 2P P x p (117) 

E c p - k = ^i 2 D+P 2 (H8) 

E* +k = ^ml + P 2 +p 2 -2Pp 3 x 3 . (119) 

To solve the restriction we introduce a further auxiliary 
variable W, which is an off-shell extension of the total 
energy |30j . by means of 

5(E; + El-E;_ q -El +k )= (120) 
J dW5{E c p + - W)S(W ~ E;_ q - E^ k ) . 

The square roots are now easier to handle two at a time, 
and the Dirac delta functions can be used to eliminate the 
two polar cosines, leaving behind only an integration over 
the auxiliary W variable and no deltas. 
The first delta 

5(E c p +E^-W) = ^6{x p -x P0 ) 

can be used to integrate over x p and fix it to 

P 2 +ml- m 2 D + W(2Ec ~ W) 
Xp " = 2~Fp 

and likewise, the second one 



(121) 



S(W - E c p _ q - E* q+k ) 



PP3 



6(x 3 - x 3o ) 



provides the x 3 integration and fixes the variable to 



P 2 + ml 



^30 



W(2E: 



p-g 



W) 



(122) 



With this kinematic work, the integrals in Eq. (|107p 
have been reduced to a four-dimensional integration 



dPdWdp 3 d(j) 3 - 



P:i 



(123) 



256^ 4 J """"pE'E^ 
x/ ff (P-p) [1 + / OT (P- P3 )] 

xJ2\M(s,t : x)\ 2 Q(k\p l ) ■ 

We employ standard Montecarlo methods in a computer 
programme to numerically calculate these integrals. In 
particular we employ the well-known VEGAS algorithm 
as coded by P. Lepage [21]. The ultraviolet integration is 
cutoff by the Bosc-Einstcin factors, and since we retain 
the pion masses there are no infrared enhancements. Con- 
vergence is rapidly achieved. 
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